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The concept of resolved virtual photons in addition to direct deep inelastic ep 
scattering is used to simulate the 2 + 1 jet - rate and the forward jet cross section, 
which cannot be described by direct LO/NLO processes. With standard DGLAP 
evolution of the parton densities of the virtual photon and the proton the HERA 
data can be described, provided a sufficiently large scale for the hard scattering 
matrix element is used. 



1 Introduction 

At HERA the measurement of the inclusive forward jet production SJl in the 
proton direction in deep inelastic events at low x and large is believed to be 
one of the ropst sensitive measurements for BFKL dynamics. The measured 
cross sectionO cannot be explained by standard deep inelastic scattering models 
and even next - to p. leading - order calculations (NLO, order a^) predict a too 
small cross sectionfl. 

y . . 2+1 

The measuredtl fractional 2 + 1 jet rate R2+1 — ^rsr, where two high pT 
jets come from the hard QCD scattering and the "+1" stands for the remnant 
jet, seems to be a factor of ~ 2 larger than predictions from standard deep 
inelastic scattering (DIS) models. 

However both measurements are rather well described by the Color - Dipole 
- Model (CDM) in its implementation in ARIADNE i. 

In this article we study to what extent the contribution of virtual resolved 
photons in addition to the standard DIS matrix element calculations (labeled 
as direct contributions) can reproduce the data. The studies arc performed 
with the RAPGAPQ Monte Carlo event generator which for the direct processes 
agrees rather well with the LEPTOB Monte Carlo program. 



2 Resolved Photons in DIS 

Resolved photon processes play an important role in photo-production, when 
high pt jets are observed. Any internal structure of the proton as well as of the 
photon can be resolved as long as the scale /.P of the hard subproccss, which is 
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of the order oip^, is larger than the hiverse radius of the proton ~ ^qcd 
and the photon l/i?^ ~ Q^. 

Resolved photon processes in DIS are calculated here using the Equiv- 
alent Photon Approximation, giving the flux of virtual transverse polarized 
photons Q'B together with the on shell matrix elements for the partonic sub- 
processes. The structure, q£ the virtual photon is given by parametrisations of 
the parton densities li3iB't3 Xjf^{xj,^i'^,Q'^), now depending on the scale 
and the photon virtuality Q^. Direct deep inelastic scattering and resolved 
virtual photon processes are simulated using the RAPGAP Monte Carlo event 
generatorl3. The hard subprocesses for resolved photons implemented in RAP- 
GAP are: gg qq, gg -> gg, qg ^ ^^r#l^ 99, QQ^ 99 and qq ^ qq. 

In DIS, resolved photon processes Ej'lj can play a role when the scale fj,^ of 
the hard subprocess is larger than the virtuality of the photon Q^: fj,^ > Q^. 
A typical NLO diagram like 7*17 — > qqg occurs also in the resolved photon 
picture, via the subprocess qg — > qg with the photon splitting 7* — > qq. When 
including resolved photon processes in a NLO calculation special care has to be 
taken to avoid double countingtJ. In leading order there is no double counting 
when adding resolved photons to direct processes since they have different final 
states (three parton and two parton). 

2.1 Scales, as and Parton Distribution Functions 

In leading order as processes the renormalisation scale iin and factorization 
scale fip are not well defined, and any choice of fij^ = fi% = , , + 
p^, and multiples thereof, is reasonable. There are essentially two competing 
effects: a larger scale gives smaller asin"^), but a larger parton density xf{x, fi^) 
at fixed x. The net effect depends on the details of the interaction. 

The effect of different choices of the scale ~ + 4 ■ p^, fi^ — + p^, 
fi^ = 4 ■ p^ , or fj,'^ — p^ is discussed in the following. 

jM^use the GRV 94 HO (DIS) parjametrisation of the proton parton densi- 
tiest3ll3 together with the SaSgam 2Dlill parametrisation for the partons in the 
virtual photon. Other choices like the GRV LO parton density of the real pltftj 
ton together with the virtual photon suppression factor of Drees - GodboleE3 
give similar results. 

3 Forward Jets 

We study the HERA data on high px "forward" jet production at low x and 
large in the context of resolved virtual photons. The cuts are taken from 
the preliminary HI analysis 0: y > 0.1, E'^ > 11 GeV, 160° < 0e < 173°, 
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Figure 1: The forward jet cross section as a function of x using the cpjts of the HI collabo- 
ration specified in the text. The data points are preUminary HI dataB. In la.) the dashed 
line is the RAPGAP prediction for direct DIS processes and the solid line is the RAPGAP 
prediction for the sum of direct and resolved contributions (using /j? = + 4p|,, GRV 
HO for the proton and SaSgam 2 DIS for the virtual photon). In 16.) the solid line is the 
RAPGAP prediction for the sum of direct and resolved contributions with fj? = 
the dashed line with fj? = ip^ and the dotted line with /i^ = p|,. The difference between 
solid and dashed lines is hardly visible. 



7° < Ojet < 20°, Ejet > 28.7 GeV, and j^?* > 3.5 GeV in the laboratory 
frame. The observed forward jet is required to have sufficiently large Ej^t and 
p:^^* in order to rely on well measured parton distributions in the proton. A 
Pt of the jet is required to be of the same order as (0.5 < p^/Q^ < 2) 
to suppress the contribution from direct LO processes which are governed by 
DGLAP evolution. The event topology is a small x deep inelastic event with an 
energetic high px jet in the forward (proton) direction. In Fig. |l|a.) preliminary 
data of the HI coUaborationB are shown together with the RAPGAPQ Monte 
Carlo prediction for standard DIS processes (direct) and the sum of direct and 
resolved processes, using the scale = + Ap^. We note that the dominant 
contribution of resolved photon processes to the forward jet cross section comes 
from qg qg. 

With a scale fi^ — + 4p|, (with pT defined in the parton parton cms) 
a perfect description of the data (Fig. |l|a^) is achieved. The effect of different 
choices of the scale /x^ is shown in Fig. |lp.). We observe no difference using 
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fi^ = Q^+p^ or — Ap^. In general a less hard scale results in a smaller cross 
section and worse agreement with the data at small x values. Using /i^ = 
the predicted cross section falls well below the measurement. This is because 
the condition /i^ = p^ > together with 0.5 < /p^ < 2 leaves only little 
space for resolved photon contributions. However with the concept of resolved 
virtual photons we obtain a improved description of the data for all choices of 
the scale /x^. 

4 2 + 1 jet rates in DIS 

Dijets (2 + 1 - jets) are searched for with a cone jet algorithm by requiring a 
p^^*" > 5 GeV in the hadronic center of mass frame, > 5 GeV^, y > 0.05, 
E' > 11 GeV and 156° < 0e < 173° according to the preliminary HI - analysis^. 
Since in that analysis p"^ > in the region < 25 GeV^, resolved photon 
interactions can be expected to play a significant role. It has been observedlj, 
that direct LO Monte Carlos like RAPGAP direct and LEPTO as weU as NLO 
calculations are unable to describe the data especially in the low region of 

< 30 GeV^ . However ARIADNE agrees with the data remarkably well. In 
Fig. ||a.) and 6.) the HI measurement U of R2+1 as a function of and x is 
shown together with the direct contribution and the sum of direct and resolved 
photons obtained with RAPGAP. The parameter setting is the same as for the 
forward jet analysis. We find remarkable agreement (using 11^ = + ^p\) 
between the data and the Monte Carlo prediction showing that especially in 
the region of phase space where < p\, resolved photons play an important 
role. Note that here also the dominant contribution from resolved photon 
processes comes from qg —>^ qg. 

We have checked R2+1 as a function of and x also for different scales 
^2 — Q2_^_p2^^ _ ^p2^ g^j-^j ^2 _ used for the resolved photon contribution. 
The contribution of resolved photons to R2+1 is larger for the largest scale 
^2 — q2 _|_ 4^p'^^ since there is a larger phase space for DGLAP evolution in 
the photon. The results are shown in Fig. and d.) Again we observe the 
same trend, that including resolved photons in DIS, the description of the 
measurements is considerably improved for all different choices of the scales 

5 Discussion 

A scale dependence is expected in all leading order calculations, but it becomes 
more important when the concept of resolved photons is applied to DIS: in 
order to resolve any internal structure of the virtual photon, the scale involved 
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Figure ^riThe 2 + 1 jet ratio -R2+1 as a function of and x. The data points are preliminary 
HI datalj. In 2a.) and 6.) the dotted line is the RAPGAP prediction for the direct contribu- 
tion, the dashed line for resolved photon contribution and the solid line is the prediction for 
the sum of direct and resolved contributions (using fi^ = + 4p|,, GRV HO for the proton 
and SaSgam 2 DIS for the virtual photon). In 2c.) and d.) the solid line is the RAPGAP 
prediction for the sum of direct and resolved contributions with fi^ = 
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has to be larger than the virtuaHty of the photon Q^: ji^ > . Using a 
large scale like — + ip^ results in a larger magnitude of the parton 
density, both on the photon and proton side, at fixed values of fractional parton 
momenta ^i, due to the larger phase space for DGLAP evolution. This effect 
is compensated to some extent by the smaller value of as at the same scale. 

In case of the forward jet cross section the stronger scale dependence be- 
comes understandable because of the requirement 0.5 < Q^/p^ < 2. Although 
here pt is measured in the lab system this cut restricts severely the phase 
space for resolved photons. 

With the concept of resolved photons, experimentally observed phenomena 
in DIS, such as the forward jet cross section and the 2+1 jet rate can be well 
described, provided a reasonable large scale for the hard scattering process 
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is used. This concept and its implementation into a full hadron level Monte 
Carlo generator can be seen as a model alternative to BFKL calculations lill. 
The next step would be to search for regions in phase space where the two 
descriptions do not agree but the data can only be described by one of them. 
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